Homopolar bond formation in ZnV 2 C>4 close to a metal-insulator transition 



00 
O 

o 



V. Pardo, 1,2 '0 S. Blanco-Canosa, 3 F. Rivadulla, 3 D.I. Khomskii, 4 D. Baldomir, 1 - 2 Hua Wu, 4 and J. Rivas 1 

1 Departamento de Fisica Aplicada, Universidad de Santiago de Compostela, E-15782 Santiago de Compostela, Spain 
2 Instituto de Investigaciones Tecnologicas, Universidad de Santiago de Compostela, E-15782, Santiago de Compostela, Spain 
Departamento de Quimica-Fisica, Universidad de Santiago de Compostela, E-15782 Santiago de Compostela, Spain 
* II. Physikalisches Institut, Universitdt zu Koln, Zulpicher Str. 77, D-50937 Koln, Germany 

Electronic structure calculations for spinel vanadate Z11V2O4 show that partial electronic delocal- 
ization in this system leads to structural instabilities. These are a consequence of the proximity to 
the itinerant-electron boundary, not being related to orbital ordering. We discuss how this mecha- 
nism naturally couples charge and lattice degrees of freedom in magnetic insulators close to such a 
crossover. For the case of Z11V2O4, this leads to the formation of V-V dimers along the [Oil] and 
[101] directions that readily accounts for the intriguing magnetic structure of ZnV^CU. 

PACS numbers: 71.27.+a;7.1.30.+h;75.25.+z 
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The transition between an antiferromagnetic (AF) in- 
sulator and a paramagnetic metal is associated to some 
of the most intriguing experimental observations in solid 
state physics. Unconventional forms of superconductiv- 
ity, electronic phase separation, and a variety of non- 
Fermi liquid behavior are well-known examples. Al- 
though not completely understood, it is however becom- 
ing clear that a strong coupling between charge, spin, 
orbital and lattice degrees of freedom is an essential in- 
gredient to explain many of these phenomena. So, it is 
important to define the properties and conditions under 
which this coupling among different degrees of freedom 
takes place. In a recent paper, Blanco-Canosa et al. [if. 
proposed that in single-valence systems, the transition 
between localized and itinerant electron behavior takes 
place through a transitional phase in which partial elec- 
tronic delocalization in the form of cation-clusters in an 
ionic matrix occurs. In this situation, strong lattice in- 
stabilities can be anticipated, due to a purely electronic 
mechanism Q. 

In this paper we explore, using ab initio calculations, 
how the effect of a small (realistic) U/W ratio naturally 
couples charge, spin and lattice degrees of freedom in a 
magnetic insulator close to the itinerant-electron limit. 
Particularly, we will discuss the case of ZnV204, whose 
magnetic and orbital properties are still a matter of de- 
bate nasi. We will show that, for small values of U, 
the most stable structure always consists of V-V dimers 
along the [Oil] and [101] directions, that helps to under- 
stand the intriguing magnetic structure of the material. 
This is a dramatic example of a strong electron-lattice 
coupling due to the partial electronic delocalization in 
the proximity of the itinerant electron limit. 

The electronic structure of the V 3+ ion (d 2 ) in a tetrag- 
onally distorted octahedral environment has an orbitally 
degenerate configuration. It is not yet clear what or- 
bital ordering, if any, could account for the experimen- 
tally found magnetic structure of the materialM] . Several 
pictures have appeared recently based on 

different considerations, that predict various possible or- 



bital orderings. An "antiferro-orbital" picture has been 
proposed Q , with the full occupation of a d xy orbital at 
each site and with an alternation of d xz and d yz orbitals 
along the c-axis. Also, a "ferro-orbital" ordering of com- 
plex orbitals (d xz ± i d yz ) [jj and an "orbital- Peierls" 
ordering^ with the orbitals in a d^-d^-d^-d^ pat- 
tern along the tetragonally compressed c-axis were sug- 
gested. Considering a fully ab initio, all-electron picture, 
including spin-orbit effects, the most stable solution is 
an alternating orbital ordering along the [Oil] and [101] 
directions of the orbitals with an unquenched orbital an- 
gular momentum d xz ± i d yz [H, such that the orbital 
angular momentum in every site is antiparallel to the spin 
moment. However, most of these works rely on the as- 
sumption that U/W is large, a purely localized picture, 
which seems to contradict recent experimental findings 
[l|. In many of these works, the lattice symmetry found 
in the existing experiments (space group I4i/amd) ficl ] 
was imposed as a rigid constraint. However, in general, 
electronic structure may have lower symmetry. If this 
were the case, it would be difficult to detect if corre- 
sponding distortions are weak, as in the case of ZnV204 
(ta g electrons involved). To take into account this pos- 
sibility, we carried out ab initio electronic structure cal- 
culations for ZnV204, valid for intermediate values of 
the Coulomb (Hubbard) repulsion U, and allowing for 
the lower symmetry, to find the optimal electronic and 
lattice structure. 

We present here full-potential, all-electron, electronic 
structure calculations based on the density functional 
theory (DFT), utilizing the APW+lo method [ll|, per- 
formed using the WIEN2k software [l2} ■ For the struc- 
ture optimization, we used the GGA (generalized gra- 
dient approximation) in the Perdew-Burke-Ernzerhoff 
(PBE) scheme [IH as an exchange-correlation functional. 
The geometry optimization was carried out minimizing 
the forces in the atoms and the total energy of the system, 
relaxing it from an initial trial state: unequal (dimerized) 
V-V distances along the [Oil] and [101] directions. For 
the electronic structure calculations we included strong 



2 




FIG. 1: (Color online) Schematic representation of the "dimer- 
ized" structure resulting from our ab initio, all-electron cal- 
culations for realistic values of U in ZnV204. Bold (thin) line 
represents short (long) in-chain bonds. The magnetic struc- 
ture is indicated by arrows. 



correlations effects by means of the LDA+U scheme [14| , 
where correlation effects are controlled by an effective U 
(U e ff= U — J), U being the on-site Coulomb repulsion 
and J the Hund's rule exchange constant. Spin-orbit 
effects have been introduced as a second variation using 
the scalar relativistic approximation [la ] . All our calcula- 
tions were fully converged with respect to the parameters 
used. 

ZnV2 04 crystallizes in a distorted cubic spinel struc- 
ture, where the V atoms form a pyrochlore lattice of 
corner-sharing tetrahedra. Because of the geometry of 
the pyrochlore lattice, the AF interactions between the 
V atoms are highly frustrated, leading to small Tn/Qcw 
ratios. 

The material undergoes a magnetic transition and 
a structural transition to a low-temperature tetragonal 
phase (c/a< 1) below 50 K [10]. The lattice formed by 
the V atoms can be described as built up by three V-V 
chains running along the [110], [Oil] and [101] directions 
(we use below the cubic setting). The magnetic struc- 
ture, found by neutron diffraction is AF along the 
[110] direction (within the ab plane), but along the [101] 
and [011] (off-plane) directions the spin ordering alter- 
nates, the V moments order TTIItt (see Fig. [1]). 

In the rest of the paper, for the sake of clarity, 
the structure derived from the tetragonal distortion de- 
scribed above will be called the "standard" structure, 
and the relaxed structure (see below) will be called the 
"dimerized" structure, to signal the formation of V-V 
dimers along the chains. 

As we said before, the anomalous variation found in the 
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FIG. 2: (Color online) Temperature dependence of the resis- 
tivity in the A 2+ V204 series. Small differences in the resistiv- 
ity for MgV2 04 and ZnV2 04 could be due to the intergrain 
boundary scattering in the ceramic pellets (their room tem- 
perature resistivity differs only by ~ 8 ilem). In any case, 
their activation energy (AE) is practically identical, as it 
should be, given the similar V-V distance in both compounds. 



pressure dependence of Tjv along the series A 2 +V 2 4 0, 
has been interpreted as the consequence of a variation of 
U/W, which first becomes progressively reduced by pres- 
sure, and then collapses close to the itinerant limit. Ac- 
cording to this, even though MgV204 and ZnV204 are 
still semiconducting, a partial electronic delocalization 
along the V-V bonds leading to the formation of short- 
range cation clusters and to a lattice instability can be 
anticipated. In Fig. [2] we present the experimental re- 
sistivity curves for different vanadates in which the di- 
valent cation at the tetrahedral site has been used to 
tune the V-V distance. The results show how the acti- 
vation energy decreases as the V-V distance does, and 
hence the metal-insulator transition is approached from 
the insulating side. On approaching this point, U/W 
is reduced progressively [l[ and hence partial delocaliza- 
tion can be anticipated along the V-V bonds. In or- 
der to check the possibility of a lattice instability due 
to this electronic effect, we performed a structural op- 
timization, using the experimentally observed magnetic 
structure (TITITI spins in [110] chains, TTIITT in [101] 
and [011] directions). As a starting point, we selected 
an artificially deformed structure that would give rise to 
V-V dimers (distances along the [101] and [011] chains 
get short-long-short-long) and let the system relax until 
forces on the atoms are smaller than 4 mRy/a.u. (small 
enough to consider the system is relaxed). The system 
relaxes to a structure away from the "standard" one, 
forming chains with an alternation of short-long V-V dis- 
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tances (Fig. [J). The same relaxed structure is obtained 
starting from different articially deformed structures, in- 
cluding the one with the chains having short-intermedite- 
long- intermediate V-V distances along the [101] and [Oil] 
chains (the structure proposed in Ref. d). This shows we 
have found a stable (its energy is lower than the "stan- 
dard" one by 30 meV/V for U e //= 0) and reproducible 
structure. 

The reason why one needs to displace the V atoms 
from their "standard" positions is because the latter is a 
local minimum for the system, as our calculations show 
(very small forces on the atoms). However, if we move 
the system slightly out of that structural local minimum, 
a different local minimum can be found, the one with 
dimerization shown in Fig. [T]with a smaller total energy 
and hence more stable. 

For carrying out this computational experiment, one 
needs to reduce the symmetry of the compound from the 
"standard" I4i/amd to the P4i2i2 space group, that al- 
lows having different V-V distances along the chains and, 
at the same time, keeping the same magnetic unit cell. 
The lattice parameters were not optimized, they were 
taken from Ref. 0- The accuracy of a GGA-based opti- 
mization of the lattice parameters is approximately 1%, 
which is worse than the experimental resolution. How- 
ever, the internal positions are much more difficult to 
obtain experimentally, specially in the case of small dis- 
tortions as the one we describe here. For that reason, the 
internal optimization we have carried out computation- 
ally becomes necessary. 

We have performed a comparison of the total energy of 
the "dimerized" structure and the "standard" one for var- 
ious values of U e ff. For this we needed to find the energy 
minimum for the "standard" structure. The V 3+ :d 2 ions 
in the tetragonal environment have an orbitally degen- 
erate configuration, with the d^j, orbital fully occupied 
and the additional electron in a d xz -d yz doublet. We 
tried different possible orbital orderings, and the energy 
minimum, for the "standard" structure, once spin-orbit 
coupling is included, is found for the second electron in 
an l z = ±1 state, with the direction of the orbital angular 
momentum antiparallel to the spin on each site, which 
agrees with previous results 0]. 

We have compared the total energy of the two struc- 
tures, the best "standard" one and the one with dimer- 
ization, for different values of U e // (see Fig. [3]). For 
small values of U e // (below about 2.2 eV) the "dimer- 
ized" structure is more stable than the "standard" one. 
Such small values of U e // would be consistent with 
our experimental findings of the system being close to 
the itinerant electron limit. It is remarkable that in 
the "dimerized" structure, even a GGA calculation is 
enough to open a tiny gap in the density of states 
(~ 0.05 eV). This occurs because in the "dimerized" 
structure, spin parallel V atoms move closer to allow 
two halves of d xz (d yz ) electrons to fill up a bonding 
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FIG. 3: (Color online)Total energies (in meV/V) of the two 
structures as a function of U e //. For small values of U e ff the 
"dimerized" structure proposed in this work is more stable. 
Note that, for any U e //, we take as zero the energy of the 
most stable structure. 



molecular orbital along the [101] ([011]) direction (see 
Fig. [1]). Then the bonding-antibonding splitting opens 
the gap and stabilizes this "dimerized" structure, which 
does not need occurrence of an orbital ordering (see be- 
low). This is contrary to the case of isostructural spinel 
MgTi204, where the structure along the off-plane chains 
is short/intermediate/long/intermediate, where no band 
gap opens within the GGA scheme [l6j |. 

For U e ff= 1.7 eV, the "dimerized" structure is more 
stable than the "standard" one by 5 meV/V. For a more 
delocalized case (smaller values of U e f / compared to the 
bandwidth W= 2 eV), the "dimerized" structure is al- 
ways the lowest-energy solution. However, if we increase 
the electron-electron on-site interaction U, i.e. driving 
the material towards the strongly localized limit, the 
most stable situation is the "standard" structure. So, 
from our calculations we conclude that as U/W is re- 
duced on approaching the itinerant behavior, a strong 
coupling between charge and lattice degrees of freedom 
takes place, determining the low-temperature properties 
of the system. The "dimerized" structure is consistent 
with the experimental evidences of lattice instabilities, 
and the fact that it is stable for small values of U e / / con- 
firms our hypotheses about its origin. It is remarkable 
that no special constraints have been introduced in the 
calculations, so the specific results obtained for ZnV204 
may be a manifestation of a general trend which may 
exist also in other systems close to a localized-itinerant 
crossover [l7j . 

In Table U we can see the V-V distances in the "dimer- 
ized" structure and we can compare them with the other 
members of the series. The shortest V-V distance is be- 
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TABLE I: V-V distances (in A) along the different directions 
in the "standard" (experimental values from Ref. 7) and the 
"dimerized" (our calculations) structure of ZnV^CU and also 
for the other members of the series AV2O4 in the cubic phase 
(our experimental results). 





in-plane 


off-plane short 


off-plane long 


"standard" 


2.98 


2.97 


2.97 


"dimerized" 


2.98 


2.92 


3.01 


CdV 2 4 


3.07 


3.07 


3.07 


MnV 2 4 


3.01 


3.01 


3.01 


MgV 2 4 


2.97 


2.97 


2.97 



low the critical distance for electron itineracy (~ 2.94 A), 
as estimated by Goodenough for a direct V 3+ -V 3+ bond 
across a shared octahedral edge in an oxide [18j . This is in 
agreement with the experimental prediction of the forma- 
tion of V-V molecular orbitals close to a metal-insulator 
transition in ZnV204 and is similar to the case of 
MgTi2C>4, where a tetramerization of the Ti chains has 
been observed [l6| . 

As mentioned above, for the "standard" structure it is 
possible to stabilize different orbital orderings. Among 
them, the most stable one has an unquenched orbital an- 
gular momentum of about 0.7 \ib per V site antiparallcl 
to the spin moment at each site. However, if we an- 
alyze the electronic structure of the "dimerized" struc- 
ture, such an orbital ordering is not found. In fact, we 
do not observe any orbital ordering. The occupations of 
the levels d xz and d yz are almost identical but the or- 
bital angular momenta are fairly small (about 0.1 ^b)- 
If we use a basis set with the real combination of orbitals: 
d xz ± dy Z , they are also equally populated. Hence, there 
is no trace of orbital ordering left once dimerization of 
the V chains occurs. The results we present here do not 
depend on a particular orbital ordering. The dimeriza- 
tion of the structure is caused by a spin-lattice coupling 
without orbital ordering being involved, the only neces- 
sary ingredient is the collapse of U/W in the vicinity of 
a metal-insulator transition. 

In all our calculations, we have assumed that the mag- 
netic structure is the one obtained experimentally 
But we have also carried out the calculations for differ- 
ent magnetic orderings to try to discern the values of the 
different exchange constants, assuming there exist an in- 
plane coupling (3i n ) and an out-of-plane coupling (3 ou t)- 
We have used the total energies obtained for various mag- 
netic couplings, with U e //= 1.7 eV, fitting them to a 
Heisenberg model (H= -JV j JijSiSj) and estimating the 
magnetic coupling. Different results are obtained in the 
"dimerized" and "standard" structures (the latter with 
an orbital ordering with l z = ±l-orbitals). In the "dimer- 
ized" structure J OM t can be subdivided in the coupling for 
the short J s and for the long bonds, J;. The in-plane cou- 
pling remains constant and is strongly AF: Ji„= -16 meV 
for both structures, mainly due to the singly occupied 



d xy orbital. Changes occur in the out-of-plane coupling: 
J out — -12 meV for the "standard" structure (highly frus- 
trated effective AF coupling). For the "dimerized" struc- 
ture, 3 S — 10 meV and J;= -3 meV. The exchange in short 
bonds is FM, and AF in the long ones. This dimerized 
structure removes the magnetic frustration and explains 
the experimentally observed off- plane TTIITT magnetic 
structure (see Fig. [I]). Our work resolves the difficulty 
in obtaining the correct magnetic structure met in other 
approaches [H, 0, 

The structure we have obtained is only stable in a limit 
close to itineracy, when U is comparable to W, and hence 
our conclusions can be applied to ZnV204 and not to 
more localized members of the series, like MnV204 and 
CdV204, because they are far from the metal-insulator 
transition. It is the closeness to the transition what leads 
to the formation of V-V dimers giving rise to the appear- 
ance of molecular orbitals with the subsequent charge 
derealization in a dimer in a material that is still semi- 
conducting. 

Summarizing, ab initio calculations show that ho- 
mopolar V-V bonding occurs in some members of the 
A 2+ V2C-4 series. The appearance of this effect is deter- 
mined by a considerable reduction of U/W, as it has been 
proposed to occur close to the itinerant-electron limit. 
Our results prove that charge and lattice degrees of free- 
dom couple strongly in magnetic insulators that approach 
the itinerant-electron limit. A possible physical picture 
explaining why there occurs dimerization in the TTIITT 
chains close to the itinerant regime is that in this case the 
exchange interaction resembles double exchange. Short- 
ening of TT bonds leads to a larger hopping and to a gain 
in double exchange energy, stabilizing this spin order- 
ing. In contrast, reduced hopping in longer V-V bonds 
weakens this tendency, allowing for TI ordering in such 
bonds. The unusual properties of many localized-electron 
systems that are close to the itinerant crossover should 
be revisited on the light of the results presented in this 
work. 

The authors thank the CESGA (Centro de Supercom- 
putacion de Galicia) for the computing facilities and the 
Ministerio de Education y Ciencia (MEC) for the finan- 
cial support through the projects MAT2006/ 10027 and 
HA2006-0119 and also the Xunta de Galicia through the 
project PXIB20919PR. F.R. also acknowledges MEC for 
support under program Ramon y Cajal. The work of 
D.Kh. and HW. was supported by DFG via the project 
SFB 608, and by the European Project COMEPHS. 



* Electronic address: vpardo@usc.es 
[1] S. Blanco-Canosa, F. Rivadulla, V. Pardo, D. Baldomir, 
J. S. Zhou, M. Garcfa-Hernandez, M. A. Lopez-Quintela, 
J. Rivas, and J. B. Goodenough, Phys. Rev. Lett. 99, 
187201 (2007). 



■5 



[2] Y. Horibe, M. Shingu, K. Kurushima, H. Ishibashi, 
N. Ikeda, K. Kato, Y. Motome, N. Furukawa, S. Mori, 
and T. Katsufuji, Phys. Rev. Lett. 96, 086406 (2006). 

[3] Y. Motome and H. Tsunetsugu, Phys. Rev. B 70, 184427 
(2004). 

[4] O. Tchernyshyov, Phys. Rev. Lett. 93, 157206 (2004). 
[5] T. Maitra and R. Valentf, Phys. Rev. Lett. 99, 126401 
(2007). 

[6] V. O. Garlca, R. Jin, D. Mandrus, B. Roessli, Q. Huang, 

M. Miller, A. J. Schultz, and S. E. Nagler, Phys. Rev. 

Lett. 100, 066404 (2008). 
[7] M. Reehuis, A. Krimmel, N. Buttgen, A. Loidl, and 

A. Prokofiev, Eur. J. Phys. B 35, 311 (2003). 
[8] D. I. Khomskii and T. Mizokawa, Phys. Rev. Lett. 94, 

156402 (2005). 

[9] S. DiMatteo, G. Jackeli, and N. B. Perkins, Phys. Rev. 
B 72, 020408(R) (2005). 
[10] Y. Ueda, N. Fujiwara, and H. Yasuoka, J. Phys. Soc. Jpn. 
66, 778 (1997). 

[11] E. Sjostedt, L. Nordstrom, and D. Singh, Solid State 



Commun. 114, 15 (2000). 
[12] K. Schwarz and P. Blaha, Comp. Mat. Sci. 28, 259 

(2003) . 

[13] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 

Lett. 77, 3865 (1996). 
[14] A. I. Liechtenstein, V. I. Anisimov, and J. Zaanen, Phys. 

Rev. B 52, R5467 (1995). 
[15] D. J. Singh, Planewaves, pseudopotentials and LAPW 

method (Kluwer Academic Publishers, 1994). 
[16] M. Schmidt, W. Ratcliff, P. G. Radaelli, K. Refson, N. M. 

Harrison, and S. W. Cheong, Phys. Rev. Lett. 92, 056402 

(2004) . 

[17] Another non-trivial possibility in this case is the sponta- 
neous charge disproportionation, although it is unlikely 
for the case of V 3+ oxides (LI. Mazin et al, Phys. Rev. 
Lett. 98, 176406 (2007)). 

[18] J. B. Goodenough, Prog. Solid State Chem. 5, 145 
(1971). 



